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Abstract 

The dispersive effects of vacuum polarization on the propagation of a strong circularly polarized 
electromagnetic wave through a cold collisional plasma are studied analytically. It is found that, 
due to the singular dielectric features of the plasma, the vacuum effects on the wave propagation in 
a plasma are qualitatively different and much larger than those in pure vacuum in the regime when 
the frequency of the propagating wave approaches the plasma frequency. A possible experimental 
setup to detect these effects in plasma is described. 
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I. INTRODUCTION 



The fast development of laser technology in the last years has enabled novel applications of 
lasers in different areas . Terawatt and petawatt table-top lasers are already employed, 
for example, to create new short- wavelength sources up to x-rays and 7-rays j^], to accelerate 
particles with enormous acceleration rates (even of the order of GeV/cm in the case of 

n n 

electrons) [4j and as drivers for nuclear fusion [5( . Strong optical laser pulses with intensities 
of order of 10 22 W/cm 2 have been already produced in laboratories [6J and intensities of 
order of 10 24 -10 26 W/cm 2 are envisaged [TJ]. In the presence of such strong electromagnetic 
fields it is becoming feasible to probe the nonlinear properties of vacuum as predicted by 
quantum electrodynamics (see the book in Ref. 
starting from the early works by Delbriick 



and the references therein). In fact, 
9], Heisenberg and Euler [kJ] and Weisskopf 
111 ], numerous theoretical papers have been devoted to investigate how the presence of 



strong; electromagnetic fields modifies the dielectric properties of the vacuum I2I, [k], Q, 



15 



161 ] . The typical electric (magnetic) field strength at which these effects are predicted 
to become apparent is given by E cr = m 2 e c'/fte = 1.3 x 10 16 V/cm (B cr = m 2 e c z /he = 
4.4 x 10 13 G) corresponding to the huge intensity: I cr = cE^/Air = 4.6 x 10 29 W/cm 2 . 
In the above expressions — e and m e are the negative charge and the mass of the electron, 
respectively. From a theoretical point of view, three different classes of electromagnetic 



fields have been considered as candidates able to "polarize" the vacuum: Coulomb fields of 
highly charged nuclei 9j, static magnetic fields 13J and laser fields 
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Q,Q, 



201 ] . So far, an experimental confirmation of the theoretical predictions has been achieved 



only for the first class with the observation of Delbriick scattering |2lJ] and the related 
process of 7-photon splitting in the field of a heavy nucleus [22J. On the other hand, the 
PVLAS (Polarizzazione del Vuoto con Laser) experiment is devoted to measure the extremely 
small ellipticity acquired by a linearly polarized probe laser after passing repeatedly through 
a vacuum region with an applied static uniform magnetic field of strength 5.5 x 10 4 G 

23J] . Experimental results have already been reported in Ref. 24j but they cannot be 
explained as a nonlinear quantum electrodynamics effect but as an effect due to the possible 
conversion of a photon into a pseudoscalar particle, called axion. Finally, in the experiments 
so far performed no interaction has been detected between strong laser beams in vacuum 

251 ] . Nevertheless, recent estimations of vacuum polarization effects (VPEs) in laser-laser 
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collisions, as in Ref. (18J in the geometrical optics limit and in Ref. [20| with the inclusion 
of diffractive effects, have shown a feasibility to observe VPEs in the near future. 

When a strong laser wave, able to polarize the vacuum, encounters matter, a plasma is 
immediately created. Already at intensities above 10 18 W/cm 2 for optical frequencies this 
plasma is relativistic and drastically modifies the laser wave evolution. Laser self-focusing, 
self- channeling and plasma self-induced transparency are only a few of the numerous effects 



arising during the propagation of a laser fie 
in Ref. 



26j . the recent reviews in Refs. 



27 



d through a relativistic plasma (see the book 
281 ] and the references therein). From another 



side, super-strong short-wavelength electromagnetic radiation fields approaching the value 
of the critical field have been predicted via the interaction of a strong laser pulse with an 
inhomogeneous plasma [29] . In this situation, VPEs could significantly modify the dynamics 
both of the plasma and of the laser field. The influence of the absorptive features of the 
vacuum polarization, namely the possibility of pair creation, on the propagation of a str ong 
electromagnetic wave through an electron-positron plasma has been considered in Ref. 30] 
where the dispersive effects of vacuum polarization have been neglected. The related back- 
reaction effects in the pair production process in crossed laser beams have been considered 
in Ref. 3l| and non-Markovian effects in the pair production kinetics have been found. 
The dispersive VPEs can significantly modify the radiation propagation in a plasma as it 



has been shown in Refs. 



32 



331 ] and more recently in Refs. 



34 



35 



36 



371 ] (see also the 



review in Ref. |28j). The authors refer to astrophysical scenarios where VPEs are due to 
a strong static magnetic field and transform the dispersion equation of a probe wave in 
such a way to allow the propagation of a new low-frequency mode in plasma. It has been 
shown also that in nonlinear regime a shock wave can be formed in magnetized plasmas 



due to the VPEs which could be important for the evolution of supernova remnants [36]. 
Also, vacuum polarization corrections to wave propagation in unmagnetized as well as in 
magnetized plasma have been investigated in Ref. 37| . Another manifestation of the VPEs 
in plasma, namely, the photon-photon scattering in a laser-produced channel in overdense 



plasma has been considered in Ref. [38( . The plasma channel traps the laser field providing 
high radiation field intensity on relatively long distances. 

In this paper we investigate the role of dispersive VPEs in the strong laser wave prop- 
agation in plasma in a new regime of physical parameters. In fact, the problem of the 
propagation of a strong laser wave in a relativistic plasma has been studied for the first time 
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in the seminal work in Ref. 



39] by Akhiezer and Polovin. In the paper in Ref. 39] an exact 



analytical solution of this problem was found in the case of the propagation of a plane wave 
with circular polarization. The solution obtained in Ref. 39J holds under ideal assumptions 
like the plane-wave one but it has provided a qualitative understanding of what happens in 
a relativistic plasma in more complicated situations. Up to now the Akhiezer-Polovin solu- 
tion remains the only exact solvable model for the strong laser wave propagation in plasma, 
while nowadays, the quantitative theoretical study of relativistic plasma dynamics in real- 
istic conditions is mostly performed via numerical solutions. In view of the ever increasing 
available laser intensities, we reconsider in the present paper the Akhiezer-Polovin problem 
by taking into account the VPEs. Since the critical intensity I cr remains many orders of 
magnitude larger than the available ones, a first-order perturbative solution of the problem 
in VPEs will be obtained. The solution shows a large enhancement of the VPEs in plasma 
with respect to those predicted in pure vacuum when the frequency of the probe field (that 
in our case is also the field that polarizes the vacuum) approaches the plasma frequency. A 
decisive role in this sense is played by the singular dielectric behaviour of the plasma near 



the plasma frequency. In contrast to this, the plasma was assumed in Refs. [32, |33|, [34|, |35( to 
be tenuous that is, the frequency of the radiation propagating in the plasma is much larger 
than the plasma frequency. In that regime the enhancement of the VPEs we consider here, 
is completely suppressed. It must also be stressed that, despite quantum electrodynamics 
is a well established theory from an experimental point of view, dispersive VPEs induced 
by laser fields have not yet been observed. The reason is essentially that the observation of 
these effects requires very intense lasers beams. In the following we propose a possible ex- 
perimental implementation of the ideas previously mentioned in order to exploit the plasma 
enhancement of VPEs to measure for the first time VPEs themselves. 

The paper is organized as follows. In the next Section the theoretical model is described. 
In Sec. IIHI exact analytical results are obtained in the case of a collisionless plasma and 
discussed with particular emphasis on the corrections introduced by the VPEs on the strong 
wave propagation. In Sec. [IV] a possible experimental setup to measure the VPEs in plasma 
is described. Finally, in Sec. [V] the effects of electron-ion collisions are taken into account 
analytically by a perturbative approach and the corresponding modifications required in the 
proposed setup are discussed. The main conclusions are summarized in Sec. IVII In the 
following, natural units with h—c — 1 will be used throughout. 
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II. THEORETICAL MODEL 



In Ref. [39| Akhiezer and Polovin solved exactly the problem of the propagation of 
a circularly polarized plane wave through a cold relativistic collisionless plasma. In their 
calculations the plasma was described as a one-component relativistic electron fluid moving 
in the presence of a fixed ion background. In the following we consider the analogous problem 
of the propagation of a strong electromagnetic field with circular polarization through a cold 



plasma but including the VPEs induced by the strong field itself. In contrast to Ref. [39| . 
we will take into account the motion of the ions because, for example, a proton moves 
relativistically in the presence of an optical laser field with intensity we are interested in of 
order of 10 23 -10 24 W/cm 2 . The approximation of cold plasma is well justified because the 
thermal motion can be neglected with respect to the highly relativistic motion driven by 
the external strong field. Finally, we also include the effects of the electron-ion collisions 
because they can play a role in the frequency region we are interested in, that is close 
to the plasma frequency. The equations that describe the previous physical scenario are 
essentially Maxwell equations in the presence of two real charged fluids with their fluids 
equations of motion, and of an additional "vacuum fluid" represented by the four-current 
(p vac (r,t), J vac (r,t)) analyzed below: 

d-E = -e(N e -ZN l )+ Pvac , (1) 

d-B = 0, (2) 

d x E + d t B = 0, (3) 

dxB-d t E = -e(iV e v e - ZNm) + J vac , (4) 

d t N e + d ■ (A e v e ) = 0, (5) 

dtNt + d • (iV>,) = 0, (6) 

dtPe + (v e • d)p e = -e(E + v e x B) - v ei m T {^ e w e - 7^), (7) 

dtPi + (vj • d)pi = Ze(E + v» x B) - Vi e m T {^i - 7 e v e ). (8) 

It is more convenient to start discussing Eqs. and (jH]). They are the continuity equa- 
tions of the electron and ion fluids having densities N\(r,t) and currents N\(r,t)v\(r,t), 
respectively: in these expressions and in the following the index A G {e, i} indicates the two 
particles species at hand, i. e. electrons and ions. The ion mass and charge are and Ze, 
respectively. Also, Eqs. ([7]) and (jSJ) are the equations of motions of the two fluids in the 



presence of the electromagnetic field E(r, t) and B(r, t) and by accounting phenomenologi- 
cally for the collisions between the two fluids. The momenta p>(r, t) are connected with the 
velocities V\(r, t) through the usual relativistic relation p A (r, t) = m A 7 A (r, t)v A (r, t) with 
7^(r,t) = (1 — t> 2 (r, t)) -1 / 2 being the relativistic Lorentz factors. The expressions of the 
relativistic dissipative terms in Eqs. (J7]) and (JSJ) can be found in Ref. 40J. In these terms 
m r = m e mi/(m e + mj) is the reduced mass and v e i and are the two effective electron-ion 
collision frequencies that will be defined and discussed in Sec. |VJ In the cold-plasma limit 
we are working in, the equations of the momenta Pa(i", t) describe completely the motion of 
the two fluids. 

The electromagnetic field evolution in the presence of the electron and ion fluids is de- 
scribed by Maxwell equations ([I])-(jl]). In these equations the VPEs are included in the 
auxiliary vacuum four-current (p vac (r, t), J vac (r, £)) appearing in Eqs. and (jlj) that is 
given by 

4a 2 

Pvac = "45^ V ' ^ - B2)E + 7(E ' B)B ] ' (9) 

Ao 2 

Jvac = {V x [2(E 2 - B 2 )B - 7(E ■ B)E] - d t [2(E 2 - B 2 )E + 7(E ■ B)B] } (10) 

where a = e 2 /4ir is the fine-structure constant. Despite the vacuum four-current is not 
made of real particles (only the electromagnetic field appears in it) it fulfills automati- 
cally the continuity equation. Maxwell equations (|T])-(j4j) with the vacuum four-current 
(Pvac( r , t), Jvac( r , £)) can be obtained by applying the variational method to the Lagrangian 
density ~% e id+«^hit where the field- matter interaction Lagrangian density J^ nt has the usual 
expression (see, e. g. Ref. 0). But, the field Lagrangian density Jzffieid is not simply the 
Maxwell one (E 2 — B 2 )/2 but it is given by 

Afield = \(E 2 - B 2 ) + ^ [{E 2 - B 2 ) 2 + 7(E ■ B) 2 ] (11) 

i. e. by the Euler-Heisenberg effective Lagrangian density that takes into account the VPEs 
at the lowest order in the limits of small field frequencies u with respect to the electron mass 

n 

and of small field amplitudes with respect to the critical fields E cr and B cr [10]. In general, 
in the presence of a plasma with density N and temperature T, the use of the effective 
Lagrangian density in Eq. (iTTj) is allowed if the frequency u is also much smaller than ksT 

m 

with ks the Boltzmann constant and N 1 ' 3 [42|. Also, in the presence of particle backgrounds 
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the Euler-Heisenberg effective Lagrangian density (TTTj) has to be corrected by addin g th e so- 



44]. 



caUed effective La^a* d enSity at «e tenure aad epical potentia, BQ 
It turns out that the finite-temperature two-loop contribution to the effective Lagrangian 
density is dominant with respect to the one-loop contribution at temperatures such that 



ksT <C fn e [45j. In fact, while the one-loop corrections are exponentially damped by a 
factor exp(—m e /kBT), the two-loop ones scale as {ksT /m e ) 4 . It can be shown that these 
terms are negligible with respect to the Euler-Heisenberg effective Lagrangian density if 
ksT -C V eE, \feB~ |46(. In the following, we will assume that all the previous conditions 
are fulfilled. Before concluding we mention that also the spontaneous electron-positron pair 
creation from vacuum has been neglected in Eqs. (JTJ)- (JHD • This approximation is very well 
justified at field amplitudes much smaller than the critical field not only at zero-temperature 
but also at finite temperatures such that ksT <C m e . This is because the one- loop thermal 



corrections to t 
shown in Ref. 



re imaginary part of the effective Lagrangian density vanish and, as it is 



451 ] for a purely electric field, the imaginary part of the two-loop thermal 



effective Lagrangian density is also suppressed by the exponential factor exp(— 7iE cr /E) like 
the imaginary part of the Euler-Heisenberg effective Lagrangian density. 



III. A SIMPLIFIED ANALYSIS OF VPES FOR THE STRONG WAVE PROPA- 
GATION IN PLASMA: COLLISIONLESS CASE 

In this Section we study how the presence of VPEs influences the propagation of a strong 
wave through the plasma. For definites, we assume that the wave electric field amplitude 
and frequency are E and lu, respectively. In general, two small parameters characterize 
our problem. One small parameter is connected with the VPEs: 5 vac = (a/457r)(E j 'E cr ) 2 . 
We remind that the effective Lagrangian density in Eq. (fTTj) is correct up to first order in 
this parameter. The second small parameter, that we indicate as 5 U , is connected with the 
electron- ion collisions. It is determined by the ratio between the collision frequencies u e i, vi e 
and the laser frequency u and it is small because the particles motion is driven by a strong 
laser field. The inclusion of collisional effects makes the analytical solution of Eqs. (CD)- (IE]) 
technically more involved. Since the conclusions of the paper will not be changed, we consider 
first the easier case of a collisionless plasma and we elucidate the physical mechanism that 
leads to the enhancement of VPEs. Then, we will discuss in detail the effects of collisions 
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among electrons and ions in Sec. |Vj 



A. Exact analytical solution of the plasma and wave equations and discussion of 
VPEs 

We study here the propagation of a transverse monochromatic electromagnetic wave along 
the positive z axis with frequency u and amplitude E and with circular polarization. If we 
set the collision frequencies u ei and v ie equal to zero in Eqs. (j7j) and (jSJ) the whole system 
of Eqs. ([I])-© can be solved exactly even including VPEs. Actually, we look for a solution 
of Eqs. ([I])-® formally identical to the Akhiezer-Polovin one: 

N e (<j ) ) = ZN i (<j ) )=N , (12) 

PeO) = -J-iM = - — (xcos0 + ysin0), (13) 
Zj to 

E(0) = £(xsin</> - y cos</>), (14) 

B(0) = nE(k cos0 + y sin 0) (15) 

where we have introduced the wave phase <fi = kz — ut = u{nz — t) with n = k/u being 
the wave refractive index. As in Akhiezer-Polovin case, we have assumed the plasma to 
remain always locally neutral. This approximation is justified for laser pulses circularly 
polarized [as that considered in Eqs. (1141) and (f!5l) ] and much longer than the plasma 
characteristic time ~ 27r/u p i t o needed to restore quasineutrality [see Eq. ([18]) below for 
the definition of the effective plasma frequency ou p i,o\- In the following we will assume that 
this is the case. In particular, we will be interested in the regime u> > u> p i t o; then we 
will implicitly consider laser pulses lasting many laser periods. The above solution has the 
property that the two electromagnetic invariants that enter the vacuum polarization four- 
current are E 2 (<p) — -B 2 (0) = (1 — n 2 )E 2 = const, and E(0) ■ B(0) = 0, respectively. In this 
way, by plugging the above expressions into Eqs. ([I])-© one easily obtains that they solve 
those equations if 




7o,i / 45tt E 2 r uj 4 \ 7o, e 70 



2 2o E . 2 \2 
n ^ + 4^El [ pl ' o) 



(16) 
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where we have introduced the non-relativistic plasma frequencies uj 2 ix = Z\N e 2 /rri\ with 
Z e = 1 and Zi = Z, the relativistic Lorentz factors 7o,a = a/1 + Z\e 2 E 2 / u 2 m? x and the wave 
refractive index in the absence of VPEs 



1 V 7o, e 7o,i /V u 



with the effective plasma frequency being 



-V.o= + (18) 

70,e 7o,i 

The above solution describes exactly the wave propagation through the plasma by including 
the VPEs. As it is clear from Eq. ffl6|) . the VPEs modify the wave refractive index. 
We observe that since we are interested in wave frequencies close to the effective plasma 
frequency uipip we are not allowed for the moment to expand further the square root in Eq. 
f fT6|) with respect to VPEs even if we have neglected from the beginning terms proportional 
to 5 2 ac with respect to those of order of unity. In this respect, it is interesting to note that, 
although we have considered the VPEs up to first order in 5 vac in the Maxwell equations 
(CE]) and (jlj), our final result contains vacuum terms to all orders because they are compared 
with quantities that in principle can be much smaller than unity. Concerning the influence 
of VPEs on the refractive index no some comments are in order. First, we observe that 
the VPEs do not result in a simple scaling of the plasma frequency. Moreover, as in pure 
vacuum, the VPEs imply an increase of the wave refractive index and then a decrease of its 
speed. One can understand intuitively this fact by noting that including the VPEs means to 
consider the possibility that during its propagation a photon can "transform" into a virtual 
electron-positron pair that then annihilates giving back the initial photon. It is sensible 
that due to this "transformation" into a pair of massive particles the photon is slowed 
down because the effective speed of the electron-positron pair is less than the speed of light. 
Finally, we point out that at small plasma densities, i. e. if n p ; o — > 1, the VPEs vanish. On 
one hand, this corresponds to the fact that in vacuum a plane wave does not polarize the 
vacuum itself |{|. On the other hand, this implies that in the situation considered here, the 
VPEs on the refractive index cannot be independent of the plasma quantities, in particular 
of its density. From this point of view we can conclude that the polarization of the vacuum 
arises here due to the essential modification of the strong wave by the electromagnetic field 
produced by the electron and ion currents in plasma driven by the strong wave itself. 



Equation f|T6l) already contains the indication that an enhancement of VPEs compared 
with those in vacuum can arise when the frequency of the laser field approaches the plasma 
frequency Thus, in the limit n p i^ — > 0, the refractive index is mainly determined by the 
VPEs: n — > y/2S va£ and it scales as the square root of the small parameter of the vacuum 
correction <5 vac . Meanwhile in vacuum, the vacuum polarization correction to the refrac- 
tive index is always proportional to <5 vac . Of course, this limiting case has little physical 
implication because of the strong reflection of the laser wave from plasma and its absorp- 
tion in plasma. Nevertheless, the mentioned tendency of enhancement of VPEs remains in 
the regime when the laser frequency is close to the plasma frequency but when the wave 
propagation in plasma still takes place. This is discussed in the next Section. 



IV. A POSSIBLE EXPERIMENTAL SETUP TO DETECT THE VPES IN 
PLASMA 



As we have seen in the previous Section, the VPEs affect the refractive index of the 
strong wave. Then, in order to detect these effects we need to make at least two waves to 
interfere. Nevertheless, due to the nonlinearity of the plasma equations, the superposition 
principle does not hold in the plasma. For this reason we consider the experimental setup 
schematically described in Fig. 1. A laser beam with intensity Iq and frequency u is 
linearly polarized along the x axis and it is split into two waves with intensities io,i and 
J 0j 2 such that I = I 01 + I 02 - Each wave passes through a quarter-wave plate and becomes 
circularly polarized but in such a way that the two electric fields rotate in opposite directions. 
Then, they enter two different plasmas with densities AT 0)1 and iVo,2- 111 the approximations 
discussed in the previous Sections we can write the electric fields of the two waves as in 
Eq. ( fl4l) by distinguishing them by means of the index j = 1,2 and by accounting for 
the fact that they rotate in opposite directions. We point out here that the electric field 
amplitude E appearing in Eqs. (|T4|) and (fT6|) is the electric field amplitude of the wave in 
the plasma. Experimentally the laser wave is produced in vacuum and then it enters the 
plasma. In the case of normal incidence the laser electric field amplitude is connected with 
the corresponding quantity in vacuum Eq by the formula E = 2Eq/(1 + no) 48( and then 
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the two electric field amplitudes Ej in the plasmas are given by 

Ej = -^—y^I^ (19) 

with no j being the refractive indices of the two plasmas. In this respect, since for each j the 
quantity no,j also depends on Ej, then both rioj and Ej must be determined consistently 
from Eqs. (fT6l) and (fT9l) once the intensity in vacuum 7 j an d the density N j are given. 
Now, if after a distance z the two waves exit the plasmas with amplitudes E' Q ■ and are made 
to interfere, the total electric field Eq(0) can be written as 

E' o (0) = Re {e«W [^(y + «) + ^, 2 (y - zx)e l2A ^>] } (20) 

where $(0) is a common phase and 2A0(;z) is the total phase difference between the two 
waves. The quantity A<j)(z) can be written as 

i n 2 _ rl 2 + 2a_]?L(T _ „2 )2 2a ®l 2 \2 
Ziffl 2 = -CJZ - . 21 

V ; 2 n 0il + n 0i2 V ; 

The wave with electric field ( 120]) is elliptically polarized but due to the phase difference in 
the second term the main axes of the ellipse are rotated with respect to the x and y axes by 
the angle A(p(z). This rotation angle is the quantity to be measured in order to detect the 
VPEs. Equation (12ip clearly shows the possibility with respect to pure vacuum of enhancing 
the VPEs in plasma by exploiting the singular behaviour of the plasma refractive index for 
laser frequencies close to the effective plasma frequency. In fact, the typical expression of 
the analogous phase difference in vacuum corresponding to 2A(p(z) is given by 2A<p vac (z) ~ 
(a/457r)(£J/ E ct ) 2 ujz. Then, if in Eq. ( 12TI) n j 1, i. e. if u > oJpipj-, the polarization 
rotation angle in plasma due to VPEs is much larger than A0 vac (z). 

Nevertheless, we have to point out that while in pure vacuum the total phase difference is 
already proportional to the VPEs, instead in plasma the phase difference contains a plasma 
contribution. Then, the problem here is to isolate the VPEs in a measurable way. To achieve 
this we require the two waves to have the same zero-order refractive index n p i t0 (N , 1 ) (for 
the sake of clarity, we have explicitly indicated the dependence of the refractive index on 
the plasma density and on the wave intensity in vacuum). Although, it is desirable that 
the laser intensities of these two beams are differing as much as possible to have a larger 
difference in the refractive index due to VPEs [see Eq. (1211) ]. We can achieve our goal by 
choosing different densities A^ 0) i and A^ 0) 2 for the two plasmas. Thus, for the fixed densities 
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No t i and Nq^ as well as the laser frequency u, the laser intensities Jo,i and Jo,2 have to be 
adjusted in such a way [see Eq. ([17]) ] that: ^,0,1(^0,1,-^0,1) = ^,0,2(^0,2,-^0,2) = n*, where 
n* is the value of the zero-order refractive index that we want to assign to both plasmas. 
Of course, in order to obtain physically "sensible" values of the intensities 7o,i an d ^0,2, 
the two densities cannot have arbitrary values. In Fig. 2 we show a numerical particular 
example of the outlined procedure. First, we have fixed the laser photon energy equal to 
u = 1.17 eV, the ion charge equal to Z = 46 (corresponding to palladium) and n* = 5 x 10~ 2 . 
We stress that with these numerical values and in the whole ranges represented in Fig. 2, 
then no,j ~ n pi,o,j — n * ■ Now, the continuous line in Fig. 2 shows the pairs (Nq,Iq) that 
solve the implicit equation n p i t0 (N , Jo) — n * '■ Instead, the dotted line corresponds to the 
approximated analytical solution [see also Eq. (TTP1) ] 

(1 + n*) 2 



Io = Io(N ) 
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e 2 N n 




(22) 



m e tf 2 (l — n* 2 ) 

obtained in the limit m« — > 00 and that, as expected, underestimates the correct intensities 
when they become of order of 10 23 -10 24 W/cm 2 at which the ion motion cannot be neglected. 
Following the above procedure, if, for example, iVo,i = 10 23 cm -3 and iVo,2 = 2 x 10 23 cm -3 
are the two plasma densities, then we would choose the two intensities as 7o,i = 7.2 x 
10 21 W/cm 2 and 7 ,2 = 3.0 x 10 22 W/cm 2 (see Fig. 2). The necessity to have high plasma 
densities comes from the fact that the effective plasma frequency u* l = uy \ — n* 2 must be 
slightly smaller than the strong laser frequency uj which lies typically in the optical domain. 
Moreover, high laser intensities decrease the effective plasma frequency. 

By means of the above described procedure, the polarization rotation angle A<p(z) gives 
a measure of the vacuum polarization effects and it is given by 

16a E 2 -E 2 (l-n* 2 ) 2 _ 2a J , 2 - J 0|1 (1 - n*f 

A(j)(z) = — — uz w — uz (23) 

457T K£ r no,i + n 0y2 457r l cr n* 

where we used the fact that, as we have said, in the previous numerical example n j w 

n pi,o,j — n * ■ By substituting these parameters in Eq. ([23]) we obtain a rotation angle 

A0(^ o ) = 6.8 x 10~ 8 rad after a propagation distance of five laser wavelengths in plasma 

corresponding to zq w 100 fim. In order to compare this result with those obtained in vacuum 

we refer to the light-by-light diffraction process described in Ref. 20|]. In fact, as we have 

observed, one single plane wave in vacuum does not give rise to VPEs. Since we have assumed 

in our example that the two beams are obtained by splitting an initial laser beam, for a fair 
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comparison we have to use in vacuum the total intensity Iq = io,i + Io,2 = 3.8 x 10 22 W/cm 2 
and we obtain A0 vac = 3.8 x 10~ 9 rad which is more than one order of magnitude less than 
A^ ). This result in vacuum can be easily derived by .ookurg at the quantity indicated 
as "0 in Fig. 2 in Ref. [20(. It shows the polarization rotation angle that an x-ray probe 
undergoes by interacting with a strong optical standing wave as a function of the distance 
between the photon polarimeter and the lasers collision region. In the figure the strong 
optical wave intensity is 10 23 W/cm 2 and, since the dependence on the polarimeter distance 
is weak and since the rotation angle is proportional to the strong optical laser intensity, the 
desired rotation polarization angle can be obtained by multiplying the average value of the 
rotation angle in the figure, which is approximately 10~ 8 rad, times Jo[W/cm 2 ]/10 23 . 

We mention here that the tuning of the intensities Iqj to eliminate the pure plasma 
contribution to A<j){z) can be avoided by performing different experiments with different 
laser peak intensities Jo- In this case, one could reveal the presence of VPEs by fitting the 
resulting experimental curve A<fi(z) = A<fi(z, Jo) by using Eq. ff2TT) with and without the 
inclusion of VPEs. 

From the previous analysis we can conclude that the presence of a plasma can give rise 
to a large enhancement of the VPEs. This result is very relevant also because, as it is 
clear from the above example, in plasma there is a significant disadvantage (that has been 
accounted for) with respect to vacuum. In fact, in vacuum one can use high-frequency, i.e. 
x-ray, probes for which the phase difference is proportionally larger than that for optical 
probes. Instead in plasma, since the probe field (that here is the same as the strong field) 
must have a frequency close to the plasma frequency one is forced to use at most optical 
probes that already require high electron densities of order of 10 23 cm -3 . 

Despite the above results, we are aware of the practical difficulties that are behind the 
experimental realization of the previous setup. The two intensities Jo,i and io,2 must be 
carefully adjusted in order to eliminate the spurious zero-order effects in A(f>(z) due to the 
presence of the plasma. Also, since the strong waves have frequencies close to the effective 
plasma frequency different kinds of plasma instabilities can arise like, for example, the 
so-called parametric resonance (see Ref. Q and the references therein). Finally, other 
difficulties have to be mentioned that are connected with the basic assumptions of our 
model. For example, the monochromatic-wave approximation has allowed us to isolate, at 
least theoretically, the VPEs and to obtain an analytical estimate of them. Nevertheless, 
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if we considered a more realistic pulsed wave with a transverse profile then the treatment 
would be much more complex and many other effects should be taken into account like laser 
self- focusing, channeling, filamentation and so on (see Ref. 27j and the references therein). 
A further remark is in order because some effects can hinder the laser penetration into 
the plasma. We mention here the parametric instabilities and the plasma density profile 
steepening. On the one hand, the excitation time of the Raman instability is extremely 



short, about half of the laser period 



491 ]. excluding the possibility to avoid it by using short 



laser pulses. Nevertheless, the parametric instabilities saturate due to nonlinear mechanisms 
before the laser beam will be depleted 50j. Moreover, the transmission of the laser wave 
through the plasma improves at laser intensities exceeding 10 20 W/cm 2 5l|], and almost 
half of the laser pulse can be transmitted through the plasma on the length of about 200 
wavelengths. On the other hand, the profile steepening of the electron and ion densities due 
to the laser ponderomotive force can play a more important role. Because of this effect the 
laser intensity threshold for induced transparency is significantly raised at a given plasma 
density, or the density threshold is decreased at a fixed laser intensity 52]. As we have 
observed, we have chosen high plasma densities of order of 10 23 cm -3 to fulfill the condition 
of the near overdense plasma uj > u p ifi with u in the optical region. Instead, we could 
exploit the density profile steepening to consider smaller initial plasma densities. In this 
case one could expect the effect of enhancement of VPEs to be present at the new resulting 
threshold of the plasma transparency but only in the steepening region. This is because the 
enhancement effect arises when the effective refractive index of the plasma tends to zero. 
In order to discuss quantitatively all the above effects a detailed analysis in more realistic 
conditions is required. This would include, for example, particle-in-cell simulations and it is 
out of the scope of the present treatment that is meant to be essentially qualitative. 

One final comment concerns the practical difficulty in measuring rotation polarization 
angles A0 with an accuracy of the order of vacuum polarization corrections, that is 10~ 8 rad 
in our numerical example. This problem is, of course, present also in the pure vacuum case 
independently of plasma. Nevertheless, very sensitive polarimeters are available today in 
the optical regime that are able in principle to measure polarization rotation angles with an 
accuracy of order of 10~ 8 rad [s3] . Also, the number of photons delivered in one pulse by 
an optical laser with intensity of order of 10 22 W/cm 2 is much larger than the number of 
photons n ~ 10 16 required from statistical considerations to measure polarization rotation 
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angles with an accuracy of order of 10 -8 rad. 

All the above analysis has been carried out by neglecting the collisional effects. This 
can be a too crude approximation mostly because we are interested in the physical regime 
where the wave frequency approaches the plasma frequency. We will include the effects of 
electron-ion collisions in the next Section and, actually, it will appear that the previous 
experimental setup has to be modified. 

V. INCLUSION OF COLLISIONAL EFFECTS 

In the previous Section we have considered the propagation of a strong wave through 
a collisionless plasma in order to point out the role of the VPEs and to show how they 
are enhanced with respect to pure vacuum and how they can be detected experimentally. 
The enhancement is significant when the laser frequency is slightly larger than the plasma 
frequency. As we have already observed, in this parameter region the inclusion of collisional 
effects is conceptually relevant. In this Section we want to investigate this case that will turn 
out technically more complex than the collisionless case already considered. Nevertheless, 
we want to point out here that the final conclusions about the enhancement and the mea- 
surability of VPEs will not be changed by the inclusion of the collisional effects. Moreover, 
the analytical implicit solution we will find in the following Paragraph [see, in particular, 
Eqs. (I45l)- (f48l) . ( 1661) and (1671) ] is of interest by itself independently of VPEs. In fact, it gives 
an accurate description of the wave propagation through the plasma in the limit of small 
collision frequencies with respect to the laser frequency. 

A. Analysis of the VPEs for the strong wave propagation through a collisional 
plasma 

Starting again from our initial Eqs. ([I])-® we can assume here that all physical quantities 
depend only on the adimensional variables tp = kz and 77 = ut. Also, we are interested in 
the propagation of a transverse electromagnetic wave. In this case it is possible, as before, 
to look for a solution of the equations such that the electron and ion charge densities are 
equal to each other and constant and uniform, and the two fluids velocities v\((p, r/) are 
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transverse: 

N e ((p,rj) = ZN i ( ( p,ri)=N , (24) 
k-v x (ip,rj) = (25) 

where k = kz is the wave vector. Since the fluids velocities are transverse and all the 
quantities depend only on time and on the space coordinate z, it follows that in Eqs. ([7]) 
and © 

(v A (^7 7 ).a)p A (^7 ? ) = 0. (26) 

Also, by multiplying Eqs. ([7j) and (jSJ) times k we conclude that v A (<£>, 77) || B(<£>, r]). Now, 
in order to guarantee the total momentum conservation during the collisions, the collision 
frequencies v e i and Ui e must satisfy the condition N e u e i = NiUi e , that is Zv e i = z/j e . From 
this relation and from Eqs. (j7]) and ([8]) it can be seen that a solution exists such that 

p e (¥^) = P(¥^). (27) 
Finally, by introducing the vector potential A(<£>, rj) such that 

E(<p,r}) = -ud v A(ip,r]), (28) 
B^^kx^.t,), (29) 

the initial set of equations (CD)- (13) reduces to the following coupled nonlinear equations for 
the quantities A(ip,rf) and p(ip,r]) (in particular, Eq. ( 1301) comes from Eq. (jlj), and Eq. 
(ED from Eq. ©): 

kx3 v B- u;d v E = -eN ( + — ) p + J vac (A), (30) 
udjjp = eudrjA — up (31) 



with 7 A = 7 A (<p, rj) = a/1 + Zlp 2 (ip,r])/m 2 x and 

rrii + Zm e 

v = v ei ■ . (32) 

rrii + m e 

For notational simplicity we have implied that the current J va c(v 9 > r ?) has to be expressed in 
terms of the vector potential A(<p, rj) [see Eqs. ffTUj) . fl28l) and fl29|) ]. 
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In order to eliminate the dependence on the variable 77, we look for a solution of the 
previous equations of the form 



A((p,7)) = A c (<f) cost] + A s (y?) sin 77, (33) 
p(PiV) = p c (y?)cos77 + p s (v9)sin77 (34) 



with the conditions 



PcO) • p,(<p) = 0, (35) 
Pc(^=PsM=P 2 M- (36) 

It is easy to see that in this way the differential equation (13 ip becomes algebraic and its 
solution is 

= 2 1 2 [^ A cM + I'A.foO], (37) 
P»(v) = 2^ 2 MM^) - ^ A c(v)]- (38) 

The momentum components expressed in this form satisfy the conditions (135]) and (|36|) if 
the analogous following conditions are satisfied by the vector potential components A c ((p) 
and A s (cp): 

A c (<p) ■ A s {cp) = 0, (39) 
A 2 M = A 2 s (ip) = A\<p). (40) 

At this point we have to solve only Eq. (130]) for the vector potential components. In fact, by 
using the above equations we can express also the relativistic Lorentz factors only in terms 



of the vector potential amplitude A(<p) as 7a (<f) = a/1 + Zl t e 2 uj 2 A 2 {ip)/m\{uj 2 + v 2 \ Now, 
the conditions (1391) and (HOI suggest to look for the vectors A c ({p) and A s ((p) of the general 
form 

A c (cp) = Ae- h ^[iccosg(<p) +ysmg(<p)), (41) 
A s (<p) = Ae- h ^[±smg(<f)-ycosg(<f)} (42) 

where the amplitude A = —E/uj is a negative quantity. With these expressions the total 
vector potential (1331) . the electron momentum ( 1341) . the electric field (1281) and the magnetic 
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field (129]) are given by 



A(<p,7i) = Ae- h Mu(<p,7i), (43) 



p F, 

pfar,) = -—e-w 



UJ 2 LOV 



(44) 



u; 2 + v 2 w ' " ' cj 2 + i/ 2 
E(^) = ^e^v^r?), (45) 
B^tj) =n J Ee-^[-/ i / (^v(^r ? )+^(^) U (^r / )] (46) 

with n = k/u. In the above expressions the prime denotes the derivative with respect to (p 
and the two perpendicular unit vectors 

u((p,7]) = ±cos(g((p) -rj) +ysin(5((v?) - ij), (47) 

v(cp, r))=x sia(g((p) -ij)-y cos(g(ip) - rj) (48) 

such that 

d v u(<p, rj) = -g'((p)v((p, rj), d v u((p, rj) = v(<p, rj), k x u(cp, rj) = -nuv((p, rj) (49) 

d^v(ip, rj) = g'(ip)u((p, rj), d v v(<f, rj) = -u(ip, rj), k x v((p, rj) = nwu(y, rj) (50) 

have been introduced. By using Eqs. f|4"5l) and f|46|) we obtain that the two quantities 

E\^-B\^=E 2 e-^\l-n 2 [h'\^)+g'\^)}}, (51) 
E(ip, rj) • B(y>, 77) = nE 2 e- 2hM ti{<p) (52) 

needed to calculate the vector J vac (^p,r]) depend only on <p. This useful property is related 
to the fact that we are considering a circularly polarized wave. It also allows, by using the 
properties of the unit vectors (|4"Tj) and (|48p . to easily eliminate the dependence on 77 in Eq. 
(130]) . For the sake of clarity we indicate here the expression of the current J va c(v 9 ) v) i n terms 
of the functions h((p) and g(<p) and of the unit vectors 11(99,77) and v((p,rj). Starting from 
Eq. (JTOl) and by using the expressions (IBTT) and (l52l) of the electromagnetic invariants, we 
obtain (for notational simplicity we omit the dependence of the various quantities on ip and 
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v) 



An 2 

Jvac = {k x d v [2(E 2 - B 2 )B - 7(E ■ B)E] - ud„ [2(E 2 - B 2 )E + 7(E ■ B)B] } 



e 

„2 



U ^- { ujEe~ h {2n 2 (E 2 - B 2 )[(h" - h' 2 + g ,2 )u + (g" - 2h'g')v] 



45mf 



+2n 2 {tiu + g'v)d v (E 2 - B 2 ) + [7nd v (E ■ B) - 2[E 2 - 5 2 )]u} 
^e-^|^e- 2ft {2n 2 [l - n 2 (h' 2 + g' 2 )][(h" - h' 2 + g' 2 )u + (g" - 2h'g')v} 

- 4n 2 (tiu + g'v) [h'[l - n 2 (h' 2 + g' 2 )} + n 2 (h'h" + #V 
+ [7n 2 (h" - 2b! 2 ) - 2[1 - n 2 (/i /2 + <?' 2 )]]u} 



(53) 



By substituting this expression of J vac (</?, 77), and the expressions (PHj) . (1451) and (SSJ) of 
p(</9, 77), E((/?,r?) and B(y9,^) into Eq. (I3"U|) and by separating the terms proportional to 
u(</5, 77) from those proportional to v(<p,rj), we obtain the following final equations for the 
unknown functions h(ip) and g(<p): 

2a E 2 



1 + 



-it \ " 2 



/ ,e^[l-n 2 (/ i ' 2 + ,' 2 )])[l + n 2 (-/ i '' + /,' 2 -,' 2 )] = ^f^ + ^ , , , 

457T E* r J ijj 1 \^e\h) li\h) J OJ A + V l 

" e- 2/l n 2 {4/i , [/i'[l - n 2 (// 2 + g' 2 )} + n 2 (/i'/i" + </(/')] - 7(h" - 2h' 2 )}, 



457T E 2 cr 



(54) 



1 + ^Z~E^ e [l-n(h +g )\\{-2hg + g) = -77T + — ttt , , 2 , „ a + 



457rE c 2 r L v ^ J ™ 2 ^ 2 \le{h) li{h)J uj 2 + u 2 

+ ^^ 2h 9'{h'[l - n\h' 2 + g' 2 )] + n\b!K> + <//)} 

(55) 

where 

/ Z 2 p 2 E 2 

7A(y) = Wl+ 2 2 ^ 2 e " 2fe(y) (56) 

are the relativistic Lorentz factors due to the particles motion in the wave field. 

In conclusion, we have reduced the initial equations ([I])-© to the above Eqs. (l54"j) and 
( 155]) . These equations cannot be solved exactly if v 7^ but we can exploit the fact that 
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in the presence of a strong wave the parameter 8 V = vjuo can be considered very sma 



give a quantitative estimate we use the usual relativistic transformation formulas 
we can write the quantity v in Eq. (152]) as 



1. To 



4l[ and 



rrii + Zm e 



mi + Zm e 



1 + 



v. 



(57) 

mi + m e mi + m e \ m e mij ,elo,i J " 

with 2/^ being the effective electron-ion collision frequency in the reference system where the 
ions are at rest. In our units v* ei is given by 54] 



7T 



3 / 2 ZN n 



4 v / 2^I(4vr) 2 m 3/ 2 ( 7 * e _ i)3/2 



log 



12tt / m\{% 



Ze 3 



Nn 



1/2- 



(58) 



where ^ = 
system with 



'1 — fQ 2 ) 1 I 2 is the electron relativistic Lorentz factor in the same reference 



'0,e 



Po 



(59) 



The index "0" in the dynamical quantities in the above expressions means that they are 
calculated by solving Eq. (f5i]) and (J55]) with v — and they coincide with those introduced 
in Sec. 1111(55]. By using the results already obtained there, it is easy to show that the 
condition 8 V <C 1 is fulfilled for electric field amplitudes E such that 



E 2 ZN oo 

> 



E 2 r ml m e 



(60) 



where for simplicity we have assumed m, — > oo and we have neglected the logarithm correc- 
tion to v. We note that the right hand side of this strong inequality is usually a quantity 
several orders of magnitude smaller than unity. For example, if Z = 10, iVo = 10 23 cm -3 
and u = 1.17 eV then ZN u/m 4 e ~ 10~ 14 . 

Now, from the zero-order solution given in Sect. II I II we realize that in the limit 5 V 1 
the functions h(ip) and g'(tp) are slowly varying because h^°'(ip) = and g^(<p) = This 
means that generally speaking h'(ip), g"(ip) ~ 5 U and h"(<p) ~ 5 2 . Moreover, as we have 
mentioned before, we have taken into account the VPEs up to first order in 5 vac already 
in the initial equations ([I])-®. For these reasons we will limit ourselves to find a solution 
of Eqs. (l54"j) and fl55l) up to first order in the two limits 5 U 1 and <5 vac < 1. In these 
approximations Eqs. (154"1) and ( |55l) assume the simplified form 



1 - nig' 2 



g - 2hg 



pl,e 



+ 



pl,i 



U 2 \le(h) 7i(/l) 
,2 



1 - 



2a E 2 e 
A5^E[ r ~ 2 



-2h 



pl,e 



+ 



U) 



pl,i 



le{h) 7i(/i) 



pl,e 



CO 



+ 



pl,i 



luo 2 V7e(^) li{h)J uo' 



(61) 
(62) 
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We observe that in these equations the two relativistic Lorentz factors are given by 
jx(h((p)) = a/1 + Z 2 e 2 E 2 e~ 2h ^ /uj 2 m\ that are different from Eq. fl56l) : for notational 
simplicity we have used the same symbol because we will not use anymore the expressions 
in Eq. ( 1561) . The substitution n = no we have tacitly made in Eqs. ( 16T|) and ( 1621) results au- 
tomatically by imposing the initial conditions h(0) = on the wave amplitude and g(0) = 
and g'{0) = 1 on the wave phase [see Eq. (16T1) ]. In turn, these conditions ensure that the 
solution of Eqs. flSTT) and (162]) reduce to the Akhiezer-Polovin solution plus VPEs if v = 0. 
In this respect, we observe that since h(0) = and h'((f ) ~ 5 U then also h(ip) ~ 5 U . Nev- 
ertheless, we have treated more exactly the terms in h(ip) in order to have a solution valid 
also for large values of (p. Now, on one hand, we obtain from Eq. (1621) that the function 
g'(<p) can be expressed in terms of the function h((f) as 



On the other hand, if we introduce the quantity 



(63) 




" (/ ' ( " n l/J ^( v70ie (^)) + 70iJ (^))J + 45vr^ ^ ( v70ie (^)) + 7o, i (^))y 

(64) 

then Eq. fIBTl) becomes 

= ^ («) 
n 



Finally, by combining Eqs. (1631) and (tool) we obtain the following differential equation for 
the function h(ip): 

1 ( U ll,e , ^pU \ 

, _ v n(h) g 7iWy , 66 x 

v ' 2uj 2 \ 7e(ft) aJ 2 rre 2 +e 2 oj 2 m 2 +Z 2 e 2 _E 2 e — 2h 

This equation can in principle be integrated and, by substituting the solution into Eq. (|65p 
we obtain 

*) = f -fv(l-=»). (W) 

In conclusion, due to the electron-ion collisions both the amplitude and the phase of the 
wave are modified. As expected, the amplitude decreases during the propagation (h(0) = 
and h'(if) > 0). Also, by neglecting for simplicity the vacuum corrections in n(h(ip)) we 
can say that if the wave frequency is larger than the non-relativistic plasma frequency 
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\Jojp le + u>p t \ then n(h((p)) > for every <p. This means that h'((f) > for every ip and then 
that the wave amplitude goes to zero for very large (p. Instead, in the complementary case 
where u> < yju^ e + ou 2 ^, the function h((p) tends asymptotically to the value h* such that 
n(h*) = that physically corresponds to the impossibility of further wave propagation. 

In order to determine the physical conditions in which the above solution holds, we 
observe that in passing from Eqs. (1541) and (155]) to Eqs. ( l6Tj) and ( 162]) we have used the 



following strong inequalities valid in the regime u > uj p i^\ 

h H {<p\ti\<p)<f/\<p\ (68) 

5 V , 5 vac < 1, (69) 

'(cp)h'(cp),5^ntg' 2 (cp)(cp)g"(cp) « S u . (70) 



Now, Eq. (16"5]) implies that g'(<p) < 1. Also, while h'(ip) ~ S v g'((p), the second derivatives 
h"(ip) and g"(<p) can be estimated from Eqs. (1661) and fl63|) as h"(ip) ~ 8 2 /n^ and g"((p) ~ 
5 u /n%. By using these estimations it can be easily shown that our approximated solution 
holds if only the two following conditions are fulfilled: 

<£«n 2 (%0), (71) 
5vac < 1. (72) 

The first of the above conditions implies that at very large (p our solution is no more valid 
because, as we have said, for wave frequencies close to u p i t o then ]im. lfi ^ 00 n(h((p)) = 0. 
Nevertheless, we can qualitatively say from Eq. ( 1661) that if we choose the initial parameters 
such that 5l <C nl then the condition (1711) is satisfied up to values of ip of order of 1/S 2 , ^> 1. 
For example, we have seen numerically that at wave intensities of order of 10 22 W/cm 2 then 
5 P ~ 10~ 10 for an optical frequency uj and, in conclusion, the condition in Eq. (ITT]) is in 
practice not restrictive at all. Finally, we have also to impose that <5 2 ac 5 V to guarantee 
that the second-order corrections with respect to VPEs are consistently negligible. This 
condition is stronger than the one in Eq. ( 1721) then our final validity-conditions are 

<5 2 «n 2 (/^)), (73) 
C « S u . (74) 

It is worth pointing out here that, independently on VPEs, the above implicit solution (166]) - 
( 167]) holds under very general conditions and it can represent by itself a useful tool to study 
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analytically the propagation of a strong wave through a plasma by including the effects of 
the electron-ion collisions. 

In order to obtain a more transparent solution we observe that if the function h(<p) itself 
is very small then we can set h((p) = in the right hand side of Eq. ([6"6"j) and h((p) ~ h'(0)(p. 
In these approximations we obtain the simplified explicit solution [see also Eq. (16 7p ] 



2 ^2 



P l,e _|_ 

V 1 70, e + 

^ ^ o„2 i 1 ( u ll,e e 2 E 2 i ^ll.i Z 2 e 2 E 2 \ 

*' v ~T 2lu 2 \ 70, e ui 2 m 2 +e 2 E 2 ^ 

70, i urni'j+Z 2 e 2 E 2 J 



Kv) - ~~ : i ..,2 ^ 7 °" ,., (75) 



9 9 x 1 / "pi, e e 2 E 2 | Z 2 e 2 E 2 \ 

\V 1 /^pi, e ^plA ^ \ 70,e a, 2 m 2 +e 2 ,E 2 ^ 70jj a,m 2 +Z 2 e 2 £ 2 J ^ 



2 a; ngu; 2 \ 7 , e 7o,« / q^ 2 i 1 f ^ e2g2 1 ^ z . 2e2£;2 

70, e aj 2 m 2 +e 2 E 2 70, < w 2 m 2 +Z 2 e 2 _E 2 



(76) 



While the approximated expression (|75l) of the function /i(<£>) holds if <C 1 instead Eq. 
(I76p requires that the more restrictive condition ip5 u <C n 2 , is fulfilled. In the following we 
will always refer to this approximated solution. To the sake of clarity we show in Fig. 3 
a typical behaviour of the amplitude of the electric field (H5|) divided by E as a function 



of the adimensional variable cp with u = 1-ly 00 pi e + u pii- The Figure clearly shows the 
exponential damping of the laser wave. The vertical line corresponds to tpa ~ 1/105„ and 
it is clear that, as predicted theoretically, the approximated solution ( i75l) works well for 



B. Modification of the experimental setup due to the inclusion of collisional effects 

Now that we have included the collisional effects in the analysis of the strong laser wave 
evolution in plasma [see Eqs. (1751) and (]76j) ]. we want to discuss which modifications are 
needed in the possible experimental setup for the detection of VPEs in plasma discussed in 
Sec. [TV] In the approximations discussed the electric field of the traveling wave is given by 
[see Eqs. (@5J and (HH]) ] 

E{z,t) = Ee- h ' {0)ujnoz [xsin (cu(n z - t) + g"(0)(con z) 2 /2) 

-ycos (u(n z-t) + g"(0)(un z) 2 /2)] . 

The refractive index uq is given in Eq. (fl6|) and the expressions of the derivatives h'(0) and 
g"(0) can be easily deduced from Eqs. (1751) and (I76|) . On one hand, the enhancement effect 
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on vacuum effect depends only on the expression of the refractive index no which, in the 
present framework, is not modified by the presence of collisional effects. On the other hand, 
the collisional effects induce a reduction of the wave amplitude and the arising of a frequency 
chirping. Actually, the attenuation of the wave amplitude does not change the conclusions 
in Sec. [TV] It is only understood that the values of the amplitudes E' 1 and E' 2 introduced 
in Eq. (1201) will be different from those used above. Instead, the total polarization rotation 
angle A$>(z) contains not only the quantity A<fi(z) already considered before [see Eq. (1211) ] 
but it also receives a chirping contribution Acp^z) that is given by 

AMz) = [#(0)nj» 2 - ^(0)ng,J (78) 

and it is proportional to z 2 . In this way, in order to isolate completely the VPEs we have to 
exclude the chirping term from the polarization rotation angle. This can be achieved, at least, 
in two ways. For example, we could consider so strong fields that |A0 c h(-z)| <C 5 vajC uz/ (n 0i i + 
^0,2)- We have shown numerically that this strong inequality is fulfilled at laser intensities 
Jo ~ 10 26 W/cm 2 for distances z of the order of ten laser wavelength in plasma. Another 
way is to perform two phase measurements A$>(za) and A<&(zb) differing only in the plasma 
lengths za and zb- In this case, from Eqs. (121]) and (ITBl) we have A<P(za)zb = A<P(zb)za 
and A4> ch (zA)z B = A(fi ch (zB)z A . Finally, these relations allow us to express the linear phase 
difference Ac/)(za) in terms of the experimental quantities A$(z^) = A<P(za) + A(p ch (z A ) and 
A$(z B ) = A(j)(z B ) + A(j) ch (z B ) as 

= A»fa)4 - ASQaK 
z B {z B -z A ) 

VI. SUMMARY AND CONCLUSIONS 

In this paper we have calculated for the first time the dispersive effects of vacuum polariza- 
tion on the propagation of a strong monochromatic wave through a cold collisional plasma. 
We have considered the easiest situation of the propagation of a circularly polarized wave 
in order to describe the new features introduced by the vacuum polarization and to obtain 
analytical estimates of the vacuum polarization effects (VPEs). The most important result 
is that, in the presence of a plasma, the VPEs are strongly enhanced with respect to those in 
pure vacuum. As we mentioned, the physical reason of this enhancement is connected with 
the singular dielectric behaviour of a plasma near the plasma frequency: the pure plasma 
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refractive index becomes much smaller than unity for a wave with frequency only slightly 
larger than the plasma frequency. We have also proposed an experimental setup to measure 
these VPEs by making two strong fields with the same frequency and different intensities 
to interfere after crossing two plasmas with different electron and ion densities. 

It is also worth noting here that, independently on VPEs, the implicit analytical solution 
we have found in Eqs. (1661) and (1671) holds under the very general conditions (1731) and (174)) 
and it can be a useful analytical tool to describe the propagation of a strong wave through 
a relativistic collisional plasma. 

We have already pointed out in Sect. [TV] the practical and conceptual difficulties of our 
approach. Nevertheless, our aim is to stress that at the high laser intensities feasible in the 
near future the inclusion of the VPEs in studying the wave propagation through a plasma 
is, on one hand, required for the sake of consistency and, moreover, it can open other and 
new possibilities to detect the VPEs themselves. 



[1] Y. I. Salamin, S. X. Hu, K. Z. Hatsagortsyan and C. H. Keitel, Phys. Rep. 427, 41 (2006). 

[2] G. A. Mourou, T. Tajima, and S. V. Bulanov, Rev. Mod. Phys. 78, 309 (2006). 

[3] P. A. Norreys, M. Zepf, S. Moustaizis et al., Phys. Rev. Lett. 76, 1832 (1996); C. Rischel, A. 
Rousse, I. Uschmann, P.-A. Albouy, J. -P. Geindre, P. Audebert, J.-C. Gauthier, E. Froster, 
J.-L. Martin, Andre Antonetti, Nature 390, 490 (1997); A. L. Cavalieri, D. M. Fritz, S. H. 
Lee et al. , Phys. Rev. Lett. 94, 114801 (2005); H. Schwoerer, B. Liesfeld, H.-P. Schlenvoigt, 
K.-U. Amthor, and R. Sauerbrey, Phys. Rev. Lett. 96, 014802 (2006); B. Dromey, M. Zepf, 

A. Gopal et al, Nature Phys. 2, 456 (2006). 

[4] S. P. D. Mangles, C. D. Murphy, Z. Najmudin et al., Nature 431, 535 (2004); C. G. R. Geddes, 
Cs. Toth, J. van Tilborg, E. Esarey, C. B. Schroeder, D. Bruhwiler, C. Nieter, J. Cary, and W. 
P. Leemans, Nature 431, 538 (2004); J. Faure, Y. Glinec, A. Pukhov, S. Kiselev, S. Gordienko, 
E. Lefebvre, J.-P. Rousseau, F. Burgy, and V. Malka, Nature 431, 541 (2004); B. M. Hegelich, 

B. J. Albright, J. Cobble, K. Flippo, S. Letzring, M. Paffett, H. Ruhl, J. Schreiber, R. K. 
Schulze, and J. C. Fernandez, Nature 439, 441 (2006); H. Schwoerer, S. Pfotenhauer, O. 
Jackel, K.-U. Amthor, B. Liesfeld, W. Ziegler, R. Sauerbrey, K. W. D. Ledingham, and T. 
Esirkepov, Nature 439, 445 (2006); J. Fuchs, P. Antici, E. d'Humieres et al., Nature Physics 



25 



2, 48 (2006); B. Hidding, K.-U. Amthor, B. Liesfeld et al, Phys. Rev. Lett. 96, 105004 (2006). 
[5] M. Tabak, J. Hammer, M. E. Glinsky, W. L. Kruer, S. C. Wilks, J. Woodworth, E. M. 
Campbell, M. D. Perry, and R. J. Mason, Phys. Plasmas 1, 1626 (1994); P. Mulser and D. 
Bauer, Laser Part. Beams 22, 5 (2004); S. Atzeni and J. Meyer-ter-Vehn, The Physics of 
Inertial Fusion (Clarendon Press, Oxford, 2004). 
[6] S.-W. Bahk, P. Rousseau, T. A. Planchon, V. Chvykov, G. Kalintchenko, A. Maksimchuk, G. 

A. Mourou, and V. Yanovsky, Opt. Lett. 29, 2837 (2004). 
[7] T. Tajima and G. Mourou, Phys. Rev. ST Accel. Beams 5, 031301 (2002). 
[8] W. Dittrich and H. Gies, Probing the Quantum Vacuum: Perturbative Effective Action Ap- 
proach in Quantum Electrodynamics and its Application (Springer, Berlin, 2000). 
[9] M. Delbriick, Z. Phys. 84, 144 (1933). 
[10] W. Heisenberg and H. Euler, Z. Phys. 98, 714 (1936). 
[11] V. Wcisskopf, K. Dan. Vidensk. Selsk. Mat. Fys. Medd. 14, 1 (1936). 
[12] J. Schwinger, Phys. Rev. 82, 664 (1951). 

[13] J. J. Klein and B. P. Nigam, Phys. Rev. 135, B1279 (1964); N. B. Narozhnyi, Sov. Phys. JETP 

28, 371 (1969); Z. Bialynicka-Birula and I. Bialynicki-Birula, Phys. Rev. D 2, 2341 (1970); 

E. Brezin and C. Itzykson, Phys. Rev. D 3, 618 (1971); I. A. Batalin and A. E. Shabad, Sov. 

Phys. JETP 33, 483 (1971); V. I. Ritus, Ann. Phys. 69, 555 (1972); J. K. Daugherty and I. 

Lerche, Phys. Rev. D 14, 340 (1976); L. F. Urrutia, Phys. Rev. D 17, 1977 (1978); J. S. Heyl 

and L. Hernquist, J. Phys. A: Math. Gen. 30, 6485 (1997). 
[14] W. Becker and H. Mitter, J. Phys. A: Math. Gen. 8, 1638 (1975). 

[15] V. N. Baier, A. I. Mil'shtein, and V. M. Strakhovenko, Sov. Phys. JETP 42, 961 (1975). 
[16] E. B. Aleksandrov, A. A. Ansel'm, and A. N. Moskalev, Sov. Phys. JETP 62, 680 (1986). 
[17] A. Di Piazza, K. Z. Hatsagortsyan, and C. H. Keitel, Phys. Rev. D 72, 085005 (2005). 
[18] T. Heinzl, B. Liesfeld, K.-U. Amthor, H. Schwoerer, R. Sauerbrey, and A. Wipf, Opt. Commun. 
267, 318 (2006). 

[19] E. Lundstrom, G. Brodin, J. Lundin, M. Marklund, R. Bingham, J. Collier, J. T. Mendonga, 
and P. Norreys, Phys. Rev. Lett. 96, 083602 (2006); E. Lundstrom, J. Collier, R. Bingham, 
J. T. Mendonga, and P. Norreys, Phys. Rev. A 74, 043821 (2006). 

[20] A. Di Piazza, K. Z. Hatsagortsyan, and C. H. Keitel, Phys. Rev. Lett. 97, 083603 (2006). 

[21] G. Jarlskog, L. Jonsson, S. Priinster, H. D. Schulz, H. J. Willutzki, and G. G. Winter, Phys. 

26 



Rev. D 8, 3813 (1973); M. Schumacher, F. Smend, and I. Borchert, Phys. Rev. C 13, 2318 
(1976); Sh. Zh. Akhmadaliev, G. Ya. Kezerashvili, S. G. Klimenko et al., Phys. Rev. C 58, 
2844 (1998). 

[22] Sh. Zh. Akhmadaliev, G. Ya. Kezerashvili, S. G. Klimenko et al, Phys. Rev. Lett. 89, 061802 
(2002). 

[23] D. Bakalov, G. Cantatore, G. Carugno et al, Nuclear Physics B (Proc. Suppl.) 35, 180 (1994). 

See also the Internet site |htt p : / / www . ts . infn . it / experiments / pvlas / pvlas .html | 
[24] E. Zavattini, G. Zavattini, G. Ruoso et al, Phys. Rev. Lett. 96, 110406 (2006). 
[25] F. Moulin, D. Bernard and F. Amiranoff, Z. Phys. C 72, 607 (1996); D. Bernard, F. Moulin, F. 

Amiranoff, A. Braun, J. P. Chambaret, G. Darpentigny, G. Grillon, S. Ranc, and F. Perrone, 

Eur. Phys. J. D 10, 141 (2000). 
[26] A. V. Borovski, A. L. Galkin, O. B. Shiryaev, and T. Auguste, Laser Physics at Relativistic 

Intensities (Springer, Berlin, 2003). 
[27] D. Umstadter, J. Phys. D: Appl. Phys. 36, R151 (2003). 
[28] M. Marklund and P. K. Shukla, Rev. Mod. Phys. 78 591 (2006). 

[29] S. V. Bulanov, T. Esirkepov, and T. Tajima, Phys. Rev. Lett. 91, 085001 (2003); N. M. 
Naumova, J. A. Nees, I. V. Sokolov, B. Hou, and G. A. Mourou, Phys. Rev. Lett. 92, 063902 

(2004) ; S. Gordienko, A. Pukhov, O. Shorokhov, and T. Baeva, Phys. Rev. Lett. 94, 103903 

(2005) . 

[30] S. S. Bulanov, A. M. Fedotov, and F. Pegoraro, JETP Lett. 80 734 (2004); Phys. Rev. E 71, 
016404 (2005). 

[31] S. Schmidt, D. Blaschke, G. Ropke, A. V. Prozorkevich, S. A. Smolyansky, and V. D. Toneev, 
Phys. Rev. D 59, 094005 (1999); J. C. R. Bloch, V. A. Mizerny, A. V. Prozorkevich, C. D. 
Roberts, S. M. Schmidt, S. A. Smolyansky, and D. V. Vinnik, Phys. Rev. D 60, 116011 (1999); 
R. Alkofer, M. B. Hecht, C. D. Roberts, S. M. Schmidt, and D. V. Vinnik, Phys. Rev. Lett. 
87, 193902 (2001); C. D. Roberts, S. M. Schmidt, and D. V. Vinnik, Phys. Rev. Lett. 89, 
153901 (2002). 

[32] Y. N. Gnedin, G. G. Pavlov, and Y. A. Shibanov, JETP Lett. 27, 305 (1978). 
[33] P. Meszaros and J. Ventura, Phys. Rev. Lett. 41, 1544 (1978). 

[34] L. Stenflo, G. Brodin, M. Marklund, and P. K. Shukla, J. Plasma Phys. 71, 709 (2005). 
[35] M. Marklund, P. K. Shukla, L. Stenflo, G. Brodin, and M. Servin, Plasma Phys. Control. 

27 



Fusion 47 L25 (2005). 
[36] M. Marklund, D. D. Tskhakaya, and P. K. Shukla, Europhys. Lett. 72 950 (2005). 
[37] J. Lundin, G. Brodin, and M. Marklund, Phys. Plasmas 13 102102 (2006). 
[38] B. Shen, M. Y. Yu, and X. Wang, Phys. Plasmas 10 4570 (2003). 
[39] A. I. Akhiezer and R. V. Polovin, Sov. Phys. JETP 3, 696 (1956). 

[40] M. Gedalin, Phys. Rev. Lett. 76, 3340 (1996). See also the book quoted therein: S. R. de 
Groot, W. A. van Leeuwen, and C. G. van Weert, Relativistic Kinetic Theory (Pergamon, 
London, 1975), pp. 18-24. 

[41] L. D. Landau and E. M. Lifshitz, Classical Theory of Fields, (Pergamon Press, Oxford, 1975). 

[42] P. Elmfors and B.-S. Skagerstam, Phys. Lett. B 348, 141 (1995). 

[43] J. I. Kapusta, Finite-temperature Field Theory, (Cambridge University Press, Cambridge, 
1989). 

[44] D. Cangemi and G. Dunne, Ann. Phys. 249, 582 (1996). 
[45] H. Gies, Phys. Rev. D 61, 085021 (2000). 

[46] The previous considerations hold at zero chemical potential. The effects of the chemical poten- 
tial are more difficult to estimate because there are no expressions available for the two-loop ef- 
fective Lagrangian at finite temperature and finite chemical potential. In this case, we proceed 
yy estimating the chemical potential as that of an ideal gas: \i = kBTln[N(2Tr /m e kBT) 3 / 2 ] 



47j. At the densities we will be interested in ~ 10 23 cm" 3 and at temperatures such that 
ksT ~ 100 <jj with lo in the optical regime [see the conditions below Eq. (jlip on the laser 
frequency] then ~ k^T and it is sensible to assume that as for the temperature also the 
effects of the chemical potential can be neglected. 
[47] L. D. Landau and E. M. Lifshitz, Statistical Physics - Part 1, (Pergamon Press, Oxford, 1980), 
p. 134. 

[48] L. D. Landau and E. M. Lifshitz, Electrodynamics of Continuous Media, (Pergamon Press, 

Oxford, 1984), p. 295. 
[49] P. Mora, Plasma Phys. Control. Fusion 43, A31 (2001). 

[50] R. K. Kirkwood, R. L. Berger, C. G. R. Geddes, J. D. Moody, B. J. MacGowan, S. H. 
Glenzer, K. G. Estabrook, and O. L. Landen, Phys. Plasmas 10, 2948 (2003); C. Rousseaux, 
L. Gremillet, M. Casanova, P. Loiseau, M. Rabec Le Gloahec, S. D. Baton, F. Amiranof, J. 
C. Adam, and A. Heron, Phys. Rev. Lett. 97, 015001 (2006). 

28 



[51] J. C. Adam, A. Heron, G. Laval, and P. Mora, Phys. Rev. Lett. 84, 3598 (2000); C. Rousseaux, 
M. Rabec Le Gloahec, S. D. Baton, F. Amiranof, J. Fuchs, L. Gremillet, J. C. Adam, A. Heron, 
and P. Mora, Phys. Plasmas 9, 4261 (2002). 

[52] F. Cattani, A. Kim, D. Anderson, and M. Lisak, Phys. Rev. E 62, 1234 (2000). 

[53] K. Muroo, M. Namikawa, and Y. Takubo, Meas. Sci. Technol. 11, 32 (2000). 

[54] H. Hora, Plasmas at High Temperature and Density: Applications and Implications of Laser- 
plasma Interaction (Springer, Berlin, 1991), p. 110. 

[55] Since we are considering also the VPEs perturbatively up to first order in <5 vac it would be 
more correct to set also e) vac = 0. Nevertheless, the VPEs do not affect the zero-order Lorentz 
factors 7 0j a that are the only dynamical quantities entering the collision frequency because 
po = m e Jjl e - 1. 



29 



Figures 



30 



A 



^ QWP 



Pol. 



Laser: 
I Q and 



I 



0,2 



TP 



0.2 



V 



FIG. 1: Schematic experimental setup described in Sec. IIV1 The "QWP" boxes represent two 
quarter-wave plates, the "iVo,i" and "iVo,2" boxes represent the two plasmas with electron densities 
iVo,i and iVo,2 and the "Pol." box the polarimeter. 
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FIG. 2: Log-Log plot of the exact numerical solution of the equation n p ifl(No, Iq) = n* (continuous 
line) and the approximated analytical solution given in Eq. (|22p (dotted line). The value of n* has 
been set n* = 5 x 10~ 2 . The two points Pi = (-/Vo,i, io,i) an d P2 = (-^0,2,-^0,2), correspond to the 
situation in the numerical example given below Eq. (|22p . 



32 




FIG. 3: Log-log plot of the typical behaviour of the amplitude of the electric field (j45]) in unit of E 
as a function of the adimensional parameter (p obtained by integrating Eq. (I66p (continuous line) 
and by using the approximated solution in Eq. (j75|) (dotted line). From Eq. (jl5j) one sees that 
|E(y?, 77) I depends actually only on ip. The laser frequency u) has been set equal to 1.1 ■J^pi e + ^ j 
. The vertical line corresponds to the value ipo ~ 1/(105^) up to which the approximated solution 

is theoretically predicted to be accurate. 
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